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Abstract 

We present a new way to compute and interpret quantum tunneling in a 1-D double-well 
potential. For large transition time we show that the quantum action functional gives an analytical 
expression for tunneling amplitudes. This has been confirmed by numerical simulations giving 
relative errors in the order of 10 -5 . In contrast to the classical potential, the quantum potential 
has a triple-well if the classical wells are deep enough. Its minima are located at the position 
of extrema of the ground state wave function. The striking feature is that a single trajectory 
with a double instanton reproduces the tunneling amplitude. This is in contrast to the standard 
instanton approach, where infinitely many instantons and anti-instatons have to be taken into 
account. The quantum action functional is valid in the deep quantum regime in contrast to the 
semi-classical regime where the standard instanton approach holds. We compare both approaches 
via numerical simulations. While the standard instanton picture describes only the transition 
between potential minima of equal depth, the quantum action may give rise to instantons also for 
asymmetric potential minima. Such case is illustrated by an example. 
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I. INTRODUCTION 



Tunneling is a characteristic feature of quantum physics, having no counter part in 

classical physics. Instantons are known to be intimately connected to tunneling. The 

physics of instantons and its relation to tunneling have been discussed in Refs. Ill [M and 

ft 

the role of instantons in QCD has been reviewed in Ref. |3|. Tunneling effects and the use 
of instantons appear in many different areas of physics, such as inflationary scenarios and 
brmation of galaxies 0, 0, Q, 7\, hot and dense nuclear matter ^|, neutrino oscillations 
fl Ej, condensed matter physics (SQUIDsl Q Q 

quantum computers based on 
superconductors Q>Q|) dynamical tunneling jly, 12 1 and chemistry (chemical bindings). 

The standard instanton picture is valid in the semi-classical regime. Infinitely many 
instantons and anti-instantons contribute to give the tunneling amplitude. In this work we 
consider the opposite regime, i.e. the deep quantum regime. We use the concept of the 
quantum action, being a kind of effective action, which takes into account quantum effects 
via tuned action parameters. The action is computed from the ground state of the system. 
In the limit of large imaginary time, the ground state wave function determines the shape of 
the quantum potential, which together with a corresponding quantum mass determines the 
quantum action. The quantum action functional then gives the exact tunneling amplitudes. 
The shape of the resulting quantum potential is different from the classical double well 
potential, i.e. it exhibits a triple well structure. We find that a double instanton (resp. 
anti-instanton) is necessary and sufficient to reproduce exactly the tunneling amplitude. 
The standard instanton approach holds for oscillator-like potentials with deep wells and 
high barriers. In contrast to that, the quantum action approach holds when the potential is 
shallow and the physics is dominated by the ground state properties (Feynman-Kac limit), 
i.e. in the deep quantum regime. In this sense the quantum action functional is a method 
complementary to the semi-classical instanton approach. 

Another approach in constructing an effective classical potential has been proposed by 
Feynman and Kleinert Though similar to the quantum action in its physical goal and 
motivation, it differs by its definition. When applied to a classical double well potential 
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|l9j it gives an effective potential different from the quantum potential. In particular, 
the quantum action generates a triple-well potential with degenerate minima (all of equal 
depth), which is not the case for the effective classical potential. 

In Sect]n]we present the tunneling model and the construction of the quantum action 
functional. Sect Mil presents numerical results on how the quantum action functional has 
been calculated and how well it fits the transition amplitude. In Sect 11 VI we compare 
numerical results from the standard instanton approach with those from the quantum 
action functional. We briefly discuss the use of the quantum action method for asymmetric 
double-well potentials in SecjV] Finally, Sect J VII gives a discussion and Sect IVIIl a summary. 



II. MODEL AND ITS QUANTUM ACTION 
A. Quantum mechanical tunneling model 

Let us consider in 1-D a classical Hamiltonian system, 

// = fUv W , (i) 

with a potential of double-well shape given by 

For simplicity, we use throughout h — c — m — 1, which makes all physicals units di- 
mensionless. The potential minima are located at ±a, a = -4=. A potential barrier of 
height B = ^ is located at x = (see Fig.JT]). The potential parameters were chosen 
such that the natural frequency of the oscillations at the bottom of each well is u — 1 for 
all A. The potential is of confinement-type, i.e. tends to infinity for large Hence the 
quantum system has a discrete bound state spectrum (no scattering states) of energies E n , 
n = 0, 1, 2, • • • . Depending on the height of the potential barrier B, the following physical 
situations may occur: (i) B < Eq. The barrier is lower than the ground state energy. 
Then the ground state wave function has a single hump at the center, (ii) E < B. The 
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FIG. 1: Classical double well potential for A = (minima at —a, +a), the ground state wave 
function and two lowest energy levels. The quantum potential V — Vo (multiplied by a scale factor 
1.7331) displays triple wells (minima at —b, 0, +b). All quantities are in dimensionless units. 

ground state energy is lower than the barrier and the ground state wave function displays 
a double hump. This case is sketched in Fig.jTj. (iii) E < E\ < ■ ■ ■ < B. The bar- 
rier is higher than the first k + 1 bound state energies. Quantum tunneling in the proper 
sense occurs in cases (ii) and (iii). The transition between regions (i) and (ii) occurs at 
A = 5.34581336- 10~ 2 . In all plots involving A, a vertical gray line represents this boundary. 
The parameter A controls the barrier height and the property of the system to be located 
in the semi-classical regime. Large A means a low barrier (quantum regime) while small A 
represents two deep (almost decoupled) wells (semi-classical regime). The tunneling am- 
plitude is given by the (imaginary time) transition matrix element corresponding to the 
transition from one minimum to the other, 

(a\exp[-HT/h]\-a) . (3) 
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B. Construction of the quantum action for large transition time 

The quantum action has been introduced in Ref. j^J. It is a local action, like the classical 
action, given by 

§[x] = [ dt ^-x 2 - V(x) , (4) 
Jo 2 

where m is called the quantum mass and V(x) the quantum potential. The quantum action 
functional is defined as a parametrisation of the quantum transition amplitude G(xji,t = 
T; Xi n , t = 0) for arbitrary transition time T in the following way 

G(xfi,t = T;x in ,t = 0) = ^ Zexp[iS[x tra j]/ft] , 

trajectories 

Xtraj ■ 5S[x tr aj] = . (5) 

Here Xtraj denotes a trajectory (stationary point) of the action S going from boundary 
point (xi n ,t = 0) to (xfi,t = T). There may be several such trajectories. 

Tunneling involves imaginary time (t — > —it). Moreover, the instanton picture of tun- 
neling is usually considered in the limit of large transition time (T — > oo). Hence let 
us consider from now on time to be imaginary and transition time T to become large 
(Feynman-Kac limit). In this limit the quantum action functional has been proven to ex- 
ist, and to give an exact parametrisation of transition amplitudes by taking into account 



only a single trajectory 



21] 



G E ucl(Xfi,t = T; X in ,t = 0) > T ^oo Z Eucl ew[-§Eud[x Eucl) /fy > ( 6 ) 

where 



T 



S Eu ci[x] = dt ^ X 2 + V(x) , (7) 

denotes the Euclidean action (following physics conventions, we dropped the overall minus 
sign occuring in the action at imaginary time. It will reappear in exp[— Seuci])- Here x^ cl 
is the trajectory which makes Se uc i stationary. We use the notation E = S[x^J[. It should 
be pointed out that, although being a stationary point of the action SeucU the trajectory 
^Euci does not necessarily always minimize the quantum action. Below we will show which 



particular trajectory must be taken in order to represent the propagator. From now on we 
drop the subscript Euclidean. 

Also in this limit, analytical relations exist between the classical potential, the quantum 
potential and wave functions (2^. For example the following relation was established j^jj, 



Here Vq denotes the minimum value of the quantum potential. Those results have been 
established on the assumption of a non-degenerate ground state and a confining potential 
with a single minimum. Tunneling involves potentials with multiple (possibly degenerate) 
minima. Therefore, we take here Eq. (jHJ) as starting point and aim to construct an expression 
which generalizes that of Eq.(jOJ). 

From the Schrodinger equation one can compute the ground state wave function ?/v( x ) 
and via Eq.fjHj) obtain the quantum potential (times the quantum mass). A plot of the 
quantum potential V(x) — Vq (up to a multiplicative factor) is shown in Fig. pp. The follow- 
ing observations can be made. First, the quantum potential, like the classical potential, is 
of confining type, i.e. it goes to infinity for large \x\. It also displays the same dependence 
in x 4 for large x. Second, there is a marked difference between the shape of the classical 
potential and the quantum potential. The former has two wells, but the latter can have 
one, two or, in the case shown in Fig. [I], three wells, located at —b, and +b. Moreover, 
Eq.fjHj) shows that V(x) = Vq, i.e. the quantum potential reaches a minimum, whenever 
the ground state wave function has derivative zero, i.e. whenever it reaches a maximum or 
a minimum. This correspondence between the extrema of the ground state wave function 
and the quantum action is a property not shared by the classical potential. Finally, the 
parity symmetry of classical potential is maintained by the quantum potential also. 

Integration of Eq. (jHj) allows to express the ground state wave function in terms of the 
quantum mass and potential, 



2m 




) 



2 



(8) 
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where Z is some integration constant. In this section, we will consider only the tunneling 
regime, where the quantum potential shows a triple well structure — other cases are simpler. 
The question arises: Which is the physically valid sign in the exponent? The answer can 
be found by looking at the shape of ip gr {x) and V(x). In the regime — oo < x < — b (regime 
I), ip gr {x) decreases when x goes from —b to — oo. In the regime — b < x < (regime II), 
ip gr (x) increases when x goes from zero to —b. In the regime < x < +b (regime III), 
il> gr {x) increases when x goes from zero to +b. In the regime +b < x < +00 (regime IV), 
i/j gr {x) decreases when x goes from +b to +00. This behavior requires the following choice 
of signs, 



ip gr (x) 



Zj exp 



Zn exp 



Zin exp 



— f b ds 



Wis) 


- Vq) 


— 00 < x < — b 


Wins) 


-Vq) 


- b < x < 




-Vo) 


< x < +b 




- Vq) 


+ b < x < 00 



(10) 



Z IV exp 

At the boundary of the regimes, the wave function has to be continuous. This implies 
relations between the factors Zi, . . . , Ziv, 



x 



-b: V4.H0 = V£H0 



gr 1 



Zj = Z n exp 



+LWf <y(*)-v°) 



x = 0: V#(0) = < 7 (0) 
x = +b: V#'(+&) = 



j 11 



Z UI exp 



+ Jo + Wf (V(s)-Vo) 



nv ■ 



Because the ground state wave function ip gr (x) and the quantum potential, Eq.(jHJ), are 
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parity symmetric, the previous equation reduces to 



^ii — tJin = 

Zi = Ziy = ZJ 

r+ b hm ~ 
J=exp + J ds\j—(V(s)-Y {] ) 



(11) 



Z is a factor which normalises the wave function to unity, and log J will turn out to be the 
quantum action of the quantum instanton. 

Taking into account those continuity conditions, Eq. (fTT|) . the wave function ip gr (x) can 
be expressed as 



ip gr (x) = Z < 



exp 



cxp 



exp 



exp 



L b(is - 



+ Jo ds 



-+b 



/Wins) - 


-V ) 


%<Y(8) - 


Vo) 


W(s) - 


V ) 




-V 



2m , 



OO < X < — b 



b < x < 



< x < +6 



+ b < x < +oo 



The terms occuring in the exponents are related to the quantum action. Because the action 
is derived from a potential, energy is conserved. In imaginary time it reads 



E = — T + V = const , 



(13) 



where f = \mx 2 tr denotes the kinetic term. Thus Eq.(|13p can be resolved for the velocity, 



i tr = ±\j^(V(x tr )-E) . 
m 



(14) 
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S x fi 



Then we can express the quantum action, 

dtf + V 

'o 

E + 2T = ET + dtmx 
'n Jo 



^ ~2 



■Si 



ET+ dx (±)J2m(V(x) - E) 



(15) 



where the sign is determined by the sign of the velocity x tr - 

Now we consider the q.m. transition amplitude in imaginary time in the limit T — ► oo 
(Feynman-Kac limit), 



G(y,T : x, 0) ~t-,oo ^r(y) exp[-£ 9r T/ft]?/y(a;) 



(16) 



The coordinates x, y may be located in any of the regimes I, II, III, IV. For example let us 
consider x,y E I (— oo < x, y < —b). Combining Eqs. (jlfipi2j) yields 

G{y,T:x,0) = Z 2 exp[-E gr T/h] 



x exp 



x exp 



ds\l^{V{s)-V Q ) 

'2fh ~ 
>I»\I^(V(8)-V ) 



y 2m ~ 
ds^{V{s)-V 0) 



o 



2m ,~ 



ds^{V{s)-V Q ) 



(17) 



There are five contributions in the exponent. Each of them can be identified with the 
quantum action of some trajectory. For example, consider the term 

—6 



exp 



It corresponds to a trajectory starting from x and approaching —b some time later, i.e., 
entering the valley of the quantum potential at —b (see Fig.jS]). Because the velocity x > 
along the trajectory, the corresponding sign in the quantum action is +. Thus the quantum 
action is 

+ / dsJ2m(V(s)-V ) . (19) 
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FIG. 2: Trajectory derived from quantum action S for boundary conditions Xj n = —2 and Xfi = 
±2 (A = i and T = 100). Full line: exp[-E]. Dashed line: exp[+S]. Minima of quantum 
potential are indicated by thin gray lines. All quantities are in dimensionless units. 

The term occuring in G is 

= exp[-is; b ] , (20) 

which is valid, when we identify E = Vq. It gives a minus sign in front of the action. 
Similarly, the other contributions to G can be expressed as some part of the trajectory 
shown in Fig.j2j. The term 

exp[-E gr T/h] (21) 

corresponds to the trajectory, where the particle rests either in the valley at —b or in the 
valley at (having identified E gr = Vq). 

To sum up, in the q.m. transition amplitude occur contributions from five different 
trajectories. They are: trajectory from x to —b, then a straight line trajectory in the valley 



cxp 



-6 



ds\l^(V(s)-V ) 
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at —5, then a trajectory going over from the valley at —b to the valley at 0, then a straight 
line trajectory in the valley 0, then a trajectory going over from the valley to the valley 
—b, then a straight-line trajectory in the valley 0, and finally the trajectory going from the 
valley —b to the final boundary point y. The transition amplitude then becomes 



G(y, T; x, 0) = Z 2 exp 



^-b 



■>s.l. 



(22) 



Note that the straight line trajectory contribution can be taken out of the action integral by 
computing the action in the potential V — Vq and adding — E gr T/H to the action afterwards. 
The full trajectory is shown in Fig. [2]. The full line corresponds the those parts of the 
trajectory corresponding to the minus sign, while the dashed line corresponds to parts with 
the plus sign. Let us take a closer look at the trajectories shown in Fig. . The trajectory 
going from x = —2 at t — to x — —2 at t — 100 has a contribution from an instanton 
and an anti-instanton. Note that this trajectory does not minimize the quantum action 
due to the contribution of the instanton anti-instanton pair. This contribution is necessary 
to correctly normalize the propagator. Other trajectories (corresponding to other initial 
and final points) are built in the same way. If we assign to each piece of trajectory a 
corresponding sign in front of the quantum action, we can write 



G(y, T; x, 0) = Z 2 exp 



^ / j ' o- x traj(v) ^traj(u) 
trajiy) 



(23) 



For other regimes, where the quantum potential has only one well, the correct trajectory 
is much simpler to find and all of it is accounted negatively in the action. 



III. NUMERICAL SIMULATIONS 



In this section we want to present numerical simulations of the quantum action and 
see how well it fits the transition amplitudes. The computation of the quantum action 
functional requires to compute a trajectory. Such trajectory is a solution of the Euler- 
Lagrange equation of motion (in imaginary time), and satisfies boundary conditions at 
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initial and final points. The numerical solution of such differential equations has been 
found to be most convenient and give stable results by using a relaxation algorithm. To 
give an example how we have proceeded, consider Fig. [2], in particular, the trajectory going 
from x = —2 to x = —2. Let us recall that we work in the regime where the transition 
time T is large. In this regime we observe that the trajectory has three pieces, where its 
motion follows the bottom of a potential valley (first from t — 5 to t — 35, second from 
t = 45 to t = 60, finally from t = 70 to t = 95). Each of those pieces of trajectory in the 
valley can be cut somewhere in the middle. As a result the whole trajectory is decomposed 
in four parts; first from t — to t — 20, second from t = 20 to t = 50, third from t — 50 
to t — 80, and fourth from t = 80 to t = 100. Each of those four pieces can be computed 
separately in a numerical way. Finally, the assignment of the sign of the quantum action 
is taken from the above theoretical analysis. 

The system provides two quantities that can be used to find a characteristic time scale, 
namely the natural frequency of oscillation in each well {T^) ale = uj^ 1 = 1) and the ground 
state energy of the system (T^ le = h/E gr ). Since E gr is in the order of 1 in the range of 
parameters considered in this paper (and since we work in units where U = 1), we simply 
take T sca i e = 1. We chose to study transitions at times T = 10 and T = 100, which are 
both large compared to the dynamical time scale. For large transition time T ^> T sca \ e 
the transition amplitude is dominated by the ground state and a few excited states. The 
wave functions and energies of those states have been computed by numerically solving the 
Schrodinger equation. The quantum action, Eq. (|4"|) . is determined by the parameter of the 
quantum mass fh and the quantum potential V(x). We computed the function m(V(x)—Vo) 
from Eq.fjHJ). In order to determine the quantum mass we need another equation. In 
the Feynman-Kac limit, one cannot directly find the quantum mass independently from 
the quantum potential. The physical reason is that in the Feynman-Kac limit physics is 
dominated by the ground state. According to Eq.Q, in the ground state wave function 
occurs the product of quantum mass and quantum potential. The mathematical reason is 
that there is a symmetry transformation which keeps the transition amplitude and 
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FIG. 3: Quantum mass vs. parameter A. T = 10. All quantities are in dimensionless units, 
the quantum action invariant. The transformation 

m — > m = m/a 

V(x) -> V'{x) = aV(x) 

T^T' = T/a (24) 

leaves the transition amplitudes G(y,T;x,0) invariant. The transformation 

rh — > m! — m/a 

V(x) -> V'{x) = aV(x) 

T^T' = T/a , (25) 

leaves the quantum action E^'J invariant. This invariance implies that the choice of the 
quantum mass m is arbitrary when T is very large: this symmetry strictly holds in the 
limit T — > oo. Numerically, we found that this freedom in the choice of m was valid to good 
precision for T > 100, that is, the value of the quantum mass had no significant influence 
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on the transition amplitudes for T = 100. However, this was found not to be the case for 
T = 10. Therefore, we need additional information about the system to find the correct 
quantum mass for T = 10. We have proceeded in the following way. For a given value of A 
we made an initial guess of the quantum mass fh. Then we computed the quantum action 
for a number of initial points X{ n and final points Xfi, taken from a set {x\, . . . ,xn}- Thus 
we generated a matrix of quantum action elements Tij. Via the quantum action functional, 
Gij = Zexp[£^], we computed a matrix of transition matrix elements. Diagonalisation of 
Gij yields eigenvalues E n T/h. Then we made a variational search in the parameter fh, until 
the energy of the first excited state Ei, obtained from diagonalisation of agreed with 
the exact value (obtained from the solution of the Schrodinger equation). On the other 
hand, previous numerical experiments have shown that the quantum mass (and also the 
parameters of the quantum potential) asymptotically converge when T — > oo. Therefore, 
we simply chose the same values of the quantum mass in the T = 100 case. The obtained 
results for the quantum mass are shown in Fig. [3] as function of the parameter A. One 
observes for large A (A > 6 • 10~ 2 ) a smooth behavior of fh. At A ~ 5.3 • 10 -2 there is a cusp. 
In our opinion the behavior for A < 5 • 10 -2 is unphysical, caused by limited numerical 
precision and the fact that in this regime the contribution to the propagator from the first 
excited state becomes non-negligeable. However, this uncertainty does not translate into 
a large error on the transition amplitudes because of their vanishing dependence on the 
quantum mass in the limit T — > oo. The relative error of the transition amplitude as a 
function of final position y, while keeping the initial position fixed (x = 0), is shown in 
Fig.|H for T = 100 and A = ^ The error is larger for large y because the numerical error 
on the quantum potential is larger for large 

IV. TUNNELING: COMPARISON WITH INSTANTON PICTURE 

We have determined the tunneling amplitude 

G Tun = G(a,T : -a, 0) = (a\exp[-HT/h] \ - a) , (26) 
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FIG. 4: Relative error E(y) 



G(y,T:0,0) QA -G(y,T:0,0) Schrod 
G(y,T:0,0) S chrod 



of transition amplitude for various 



final positions y for T = 100 and A = All quantities are in dimensionless units. 

where ±a are the minima of the classical potential, for various values of the potential 
parameter A and transition time T. As reference value, the tunneling amplitude has been 
computed by solving the Schrodinger equation. Next, the tunneling transition amplitude 
has been computed from the quantum action functional, that is by computing the quantum 
potential and the quantum mass for each of the different classical potentials and transition 
times and by using the tools developed in SectlH] Moreover, the transition amplitudes have 
been obtained by using the semi-classical multi-instanton expression at two-loop order, 
given by Ref. [3|], 



G 



Tun 



8Sc. 



exp 



uoT 



1 



3S n 



sinh 



'6Sn 



7T 



exp 



— So 



71 J_ 

725^ 



(27) 



where K — 1 and So = 1/(12A) is the action of the instanton between classical minima. 
Results are shown in Figs. |5I6] for T = 10 and T = 100, respectively. The relative differ- 
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FIG. 5: Tunneling amplitude G(a, T; —a, 0) vs. potential parameter A. T = 10. Comparison 
of exact result (from Schrodinger eq.) with result from quantum action functional and from 
instanton method (semi-classical). All quantities are in dimensionless units. 

ence between the exact tunneling amplitude (from Schrodinger equation) and the quantum 
action amplitude for T = 10 and T = 100 is shown in Fig. [7]. The instanton method 
is valid in the semiclassical limit where the action So of the classical instanton is large 
compared to unity (corresponding to small A). This is clearly the case in Figs |5l6j where 
the semi-classical tunneling amplitude is quite close to the exact amplitude for small A. In 
contrast, the tunneling amplitude computed from the quantum action is better at large A. 
It diverges substantially from the exact result for small A. This divergence is not due to 
a bad choice of quantum mass. It can rather be explained by the fact that the quantum 
action is constructed from the exact propagator in the Feynman-Kac limit, where only the 
ground state contributes to the physics of the system. However, when A decreases, the 
wells get deeper and the energy of the first excited approaches the ground state energy. 
Then the first excited state contribution cannot be neglected anymore (the ground state 
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FIG. 6: Same as Fig.jS] but T = 100. All quantities are in dimensionless units. 

becomes degenerate). We can verify that the quantum action functional indeed reproduces 
correctly the ground state contribution to the propagator \ip gr )e~ EgrT ^ h (ipgr\ fc> r ah A, up 
to numerical limits. The relative error of tunneling amplitudes from the quantum action 
is shown in Fig. [Jj. The fact that the error is generally higher for T = 10 confirms this 
observation (lower transition time means that the first excited state is more important in 
the transition amplitude for the same value of A). It turns out that the quantum action 
can be used to compute transition amplitudes for a wide range of potential parameters. 

V. TUNNELING IN ASYMMETRIC POTENTIALS 

Another advantage of the quantum action functional is that it can be used to study 
tunneling in an asymmetric double-well potential, with wells of different depth. While 
such a system is much harder to study in the standard instanton picture, it does not 
require substantially more work in the quantum action context. In fact, the quantum 
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FIG. 7: Relative error of the tunneling amplitude from the quantum action functional with respect 
to exact result, for T = 10 and T = 100. All quantities are in dimensionless units. 



potential will still display one, two or three minima, having all the same depth. Instantons 
trajectories similar to those in the standard double-well case then appear likewise. For 
example, consider the potential 

2 



V(X) = 50 



X 



+ 



1 

250 



5 

x 

2 



(28) 



which is shown in Fig. jS] . This potential has two wells, but the left well is higher than the 
right one, and therefore a classical instanton does not exist. The quantum potential (up 
to a multiplicative factor) recovered from the ground state wave function is also shown in 
Fig.jSj. We find that the quantum mass in this case (for T = 100) is about rh = 0.11 and 
that the tunneling transition amplitude G(a, T; —a, 0) reproduces the exact value with a 
relative error of 2.87 • 10 -4 . 
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FIG. 8: Asymmetric double-well potential, Ea. (j28|) . the ground state wave function and corre- 
sponding quantum potential V — Vo (up to a multiplicative factor). All quantities are in dimen- 
sionless units. 

VI. DISCUSSION 

Using the quantum action functional to describe tunneling in a double well potential 
gives good results in the deep quantum regime (opposite to the semi-classical regime). 
Quantum effects appear in the quantum action via tuned parameters. In the case of tun- 
neling in a classical double-well potential this leads to a quantum potential with a different 
well structure, having possibly one, two or three wells. Like in the standard instanton pic- 
ture, quantum instantons occur also in the tunneling amplitude obtained by the quantum 
action functional. Instantons play an important role in many domains of physics so one 
may ask how important are the instantons in this case. First, we saw that they play a 
normalization role in the trajectories, that is, they ensure that the propagator is correctly 
normalized. Second, as much as standard instantons represent the only finite action solu- 
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tion (in the limit T — > oo) contributing to the propagator between the classical minima, we 
can interpret the new instantons as the only finite action solutions (in the T —>■ oo limit) 
contributing to the propagator between the minima of the quantum potential. Because these 
minima of the quantum potential correspond in position to the extrema of the ground state 
wave function, the double quantum instanton going from —b to to b represents the trajec- 
tory reproducing the largest transition amplitude. In particular, the tunneling amplitude 
between the minima of the quantum potential (— b\ exp[— HT/h]\b) is larger than the tun- 
neling amplitude between the minima of the classical potential (— a\ exp[—HT/ h]\a). Third, 
one might wonder why in the quantum action functional only two instantons, respectively 
anti-instantons, contribute, while in the standard approach infinitely many instantons and 
anti-instantons are needed. The answer is simply that the quantum action was explicitely 
constructed to reproduce exactly, in the Feynman-Kac limit, the transition amplitude, using 
a single trajectory. This is sufficient to discard the use of multiple instanton trajectories. 
Fourth, let us compare the structure of the classical instanton with the quantum instanton. 
The instanton solution corresponding to the classical Hamiltonian of Eqs. ()ll2|) is given by 



x inst {t) = — = tgh 
v 8A 



--(t-t c ) 



f&B tgh 



'Am 



(t - U) 



(29) 



'4m 

where B denotes the barrier height of the potential and t c is the center of the instanton 
The steepness of the instanton at its center is given by 



d , s 2B 

-TL x inst\t = t c ) = \ • (30) 

at V m 

Fig.p shows that that the barrier height of the quantum potential is much smaller than 
that of the classical potential. The quantum mass is of the same order of magnitude as 
the classical mass. From Eq. (j30|) we expect that the steepness of the quantum instanton is 
smaller than that of the classical instanton. This is confirmed by a numerical calculation, 
comparing the classical with the quantum instanton, shown in Fig.jU]. Also because the 
location of the potential minima are closer for the quantum potential than for the classical 
potential, the quantum instanton has a smaller action than the classical one. The bottom 
line is that the quantum instanton is "softer" than the classical instanton. It is interesting 
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FIG. 9: Comparison between the classical and quantum instantons, for A = A. All quantities are 
in dimensionless units. 

to note that a similar observation has been made previously in the context of comparing 
classical chaos with quantum chaos, again using the quantum action functional [23]. It has 
been found that the quantum action yields a less chaotic phase space than the classical 
action. The underlying reason for such behavior is unknown to us. We believe that a 
promising strategy may be to analyze the path integral and its relation to the quantum 
action functional. 

VII. SUMMARY 

This work is about tunneling described in terms of the quantum action functional. This 
point of view is complementary to the standard instanton picture: While the latter holds in 
the semi-classical regime, the former holds in the deep quantum regime. The observation 
that the quantum potential has new minima (in number and location) beyond those of the 
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classical potential may be of interest for cosmology and inflationary models. There are 



experiments on tunneling in condensed matter, e.g. Josephson junctions, SQUIDS 



12 



3 

or in atomic physics in dynamical tunneling of atoms in a time- dependent exterior field 
[lrl \v\ . It would be interesting to explore if the quantum action functional can be applied 
to describe such physics. This requires further development, in particular, to explore the 
quantum action functional in real time and for explicitely time-dependent systems. Likewise 
one may ask if tunneling out of meta-stable states can be described in our approach. We 
have shown that an asymmetric double-well potential which may give rise to meta-stable 
states (which occur, e.g., in nuclear fission and emission of a-particles) can be treated 
in this framework, provided that the potential is bounded from below (which excludes a 
potential like V ~ x 3 ). However, the description of tunneling from quasi- stationary states 
to the ground state or other excited states would require to apply the quantum action in 
real time. This will be a subject of further studies. 
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